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Abstract 

Boundary conditions are derived that determine the penetration of 
spin current through an interface of two non-collinear ferromagnets with 
an arbitrary angle between their magnetization vectors. We start from 
the well-known transformation properties of an electron spin wave func- 
tions under the rotation of a quantization axis. It allows directly find 
the connection between partial electric current densities for different spin 
subbands of the ferromagnets. No spin scattering is assumed in the near 
interface region, so that spin conservation takes place when electron in- 
tersects the boundary. The continuity conditions are found for partial 
chemical potential differences in the situation. Spatial distribution of 
nonequilibrium electron magnetizations is calculated under the spin cur- 
rent flowing through a contact of two semi-infinite ferromagnets. The 
distribution describes the spin accumulation effect by current and corre- 
sponding shift of the potential drop at the interface. These effects appear 
strongly dependent on the relation between spin contact resistances at the 
interface. 



1 Introduction 

A branch of the solid state physics and electronics called "spintronics" U devel- 
oped rapidly last years. The name is due to the decisive role which the electron 
spin and related magnetic moment play in the transport phenomena studied. 
The most important phenomena appear in magnetic junctions containing fer- 
romagnetic layers. Disturbance of spin equilibrium occurs when spin-polarized 
current flows through the interfaces of the layers. The disturbance leads to a 
number of new spin dependent contact phenomena, which are of interest for a 
theory development and for using under interpretation of experimental data. 
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To solve equations of motion for the junctions, we should have true boundary 
conditions at the interfaces. The problem of boundary conditions appears here 
by a natural way. The most of the published works consider either a contact 
between ferromagnetic and nonmagnetic materials or between two ferro- 

magnets with collinear magnetic moments such as a domain wall |SJ ■ In the both 
cases, there is a single quantization axis. Meanwhile, the boundary conditions 
determining spin current through a contact of two noncollinear ferromagnets are 
significant for a number of problems concerning spin-polarized current induced 
spin switching 7 ■ Such type of boundary conditions are treated in the present 
work. 

Using the boundary conditions derived, we calculated further the following 
contact phenomena: magnetization distribution, spin accumulation shift of the 
contact potential drop and spin accumulation contribution to contact resistance. 



2 Boundary conditions for magnetization flux 

The electron magnetization distribution rrii{x, t) (i = x, y, z) is described by 
the continuity equation 

drrii rrii-fhi 

-7T— + Vfe Jik = , (1) 

ot T 

where fhi is equilibrium value of the electron magnetization, r is spin relaxation 
time, Jik is the electron magnetization flux density (the first index determines 
the magnetization vector direction in the flux, the second one indicates the flux 
propagation direction). 

Basing on the well known derivation of the quantum mechanical formula for 
the particle (electron) current density |S] , we obtain the electron magnetization 
flux density in the following form: 

P, S1,S 2 

-^MVk^M), (2) 

where ^> PjS (r, t) is electron wavefunction with momentum p and spin state s, 
[iB is Bohr magneton, er = {a x , a y , <j z } are Pauli matrices. 

Let us see how the magnetization flux transforms under rotation of the quan- 
tization axis. Such a rotation can be due to electron transfer from one to another 
magnetic layer of the junction with different direction of the quantization axis 
as well as rotation of a homogeneous medium quantization axis (by an external 
magnetic field, for example). 

Let electric current flow along x axis with quantization axis parallel to z axis. 
The magnetization flux density has single component J zx . Then the electron 
current goes into another layer with quantization axis parallel to z' axis that 
makes an angle \ with z axis. The quantization axis rotation is described by a 
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spin wave function transformation matrix [8] 



fj r _A _ ( cos(x/2) isin(x/2)\ 
U *W- [ism( X /2) cos( x /2) ) 



(3) 



Such a transformation of the wave functions leads to transformation of the 
magnetization flux density, so that a longitudinal component J z > x — J zx cos x 
appears with polarization along the new quantization axis z' as well as a trans- 
verse component J V ' X — J zx sinx with perpendicular polarization. Very differ- 
ent spin relaxation times correspond to the longitudinal and transverse polar- 
izations, so only the longitudinal component J z > x — J zx cos x survives beyond 
a thin layer of Fermi wavelength thickness (the so called Berger-Slonczewski 
layer, see [7] for details). This gives a boundary condition for the electron 
magnetization flux density at the interface between the ferromagnets x = 0: 



Jzx cos x 



J z' X 

x=-0 



(4) 

.T = + 



Let us see how the partial spin-polarized current densities for spin-up and 
spin-down electrons transform under changing the quantization axis. We have 

3+ + 3- = 3, (5) 
e 

3+-3- = — Jzx, (6) 

Mb 

where j is total current density. From Eqs. |J5J and 10 we obtain 

]± = 1 -( 3 ±—Jzx]. (7) 



2 V Ms 

The current density j does not change under the quantization axis rotation, 
while the magnetization flux density J zx transforms in accordance with Eq. . 
Therefore, the transformed partial current densities take the form 

4 = i(.±^cos X ). (8) 

With Eqs. J5J and @ taking into account, a transformation law for the partial 
current densities takes the form 

3± = 3±coa -+j T am -. (9) 

The electric current transformations © were obtained previously in P~[T^] 
by other ways. 



3 Boundary conditions for chemical potentials 

The other boundary condition is imposed on the partial chemical potentials 
£± corresponding to spin-up and spin-down electrons. In its derivation, we 
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start from the energy flux continuity condition at the interface between two 
ferromagnets. 

The Boltzmann equation for spin-up and spin-down electron distribution 
functions / p ± takes the form 



df P ± . , ,df P ± dip(r)df p ± 

+ v±(p)^£ r -e— £-=J p±) 



Ot 



Or 



dr dp 



(10) 



where <p(r) is electric potential, v±(p) = de p ±/dp is electron velocity, e p ± is 
the electron energy with momentum p, I p ± is collision integral. 

Multiplying Eq. (|10f> by electron energy e p ± and summing over p gives 
energy conservation law 



dU± 
dt 



divWj 



P ± i p±, 



(ii) 



where U± = J2 £ p±fp± are partial electron energy densities, 



w± = ( £ P ± + e( p) v ± (p) / P ± = -j± (c± + e( p) 



(12) 



are partial energy flux densities for completely degenerate electrons. 

Consider a contact in x — plane between two homogeneous semi-infinite 
ferromagnets with different quantization axes. The total energy flux density 
along x axis is 



W = W + + W. 
1 



1 



[]+ {(+ + eip)+u (C-+ep)] 



9 3 (C+ + C- + 2e V ) + — 
2e 2/ib 



Jzx (C+-C-), 



(13) 



where J zx is the electron magnetization flux density; we used Eq. Q here. 
With boundary condition Q taking into account, the energy flux continuity 



condition W\ 



= W 



x=+Q 



takes the form 



j|(C++C-+2e^) 
— J Z J ((+-(- 

MB I 



= _ Q -(C + +C+2e V >) 



x=+0 



x=+0 ^ 

cos x \ = 0. 



(14) 



The form of the boundary conditions should not depend on values of the 
current and magnetization flux. Therefore, the contents of both curly brackets 
in Eq. I|14[) are to vanish each separately. This gives the following boundary 
condition for the partial chemical potentials at the interface: 



(C+-C-) 



x=-0 



(c; - c) 



x=+0 



cos x- 



(15) 
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From Eq. I|15|). a boundary condition can be obtained for nonequilibrium 
electron magnetization Am = (1b (An + — An_), where An± are deviations of 
the partial electron densities in spin subbands from their equilibrium values 
n±. Because of the neutrality condition we have An + + An_ = 0, so that 
An± = ±Am/2/ifl. The partial chemical potentials are related with An±, 
namely, 

- An± Am - An' ± Am' 

C±-C= = ±7. , C±-C=— — = ±~ r> ( 16 ) 

g± 2fj, B g± g± ^b9± 

where g±, g' ± are the densities of states at the Fermi level for spin-up and 
spin-down electrons, (, C are the equilibrium chemical potentials of two ferro- 
magnets. 

With Eq. Ijl6|l taking into account, we find 

C+-C-=NAm, N = -?-(— + —) . (17) 

2/"b \g+ g-J 

From Eqs. i|15fl and (|17() . we obtain the following boundary condition for 
nonequilibrium electron magnetization: 

NAm\ x= _=N'Am'\ x=+Q cos X - (18) 

The boundary conditions 10}, (|15|) and (|18|) correspond to current flowing in 
the positive direction of x axis. Under opposite current direction, the conditions 
may be analogously presented in the form 

Jz X \ x= _ = J z 'x\ x=+0 cos X , (19) 

(c + -c-)L = _ co SX = (c;-C)L =+0 , (20) 

TV Ami n cosv = N'Am'\ ... (21) 

la; — — ^ la;— +0 v ' 

4 Electron magnetization distribution 

As an illustration, let us consider spin flux transfer through a contact between 
two semi-infinite noncollinear ferromagnets in plane x = 0. The electron mag- 
netization distribution is described by Eq. with magnetization flux density 

J„ = !±Qj (22) 

e ox 

where D = {D + tj- +Z?_er + )/ a is effective spin diffusion constant, D± are partial 
electron diffusion constants, a± are partial conductivities, a = o~+ + c_ is total 
conductivity, Q = (o~+ — ct_ ) j a is conduction polarization coefficient (see 
for details). The stationary solution of Eq. with boundary conditions 
and (|18fl takes the form 

Am(x < 0) = /i_b?i ^ exp(£c/t), (23) 

3d v + cos z x 
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Figure 1 : The nonequilibrium magnetization distribution (in dimensionless vari- 
ables) at different values of the spin resistance ratio v. 1 — v = 10 (red), 2 — 
v = 1 (green), 3 — v = 0.1 (blue). Q = 0.35, Q' = 0.15, 1 = 1', n = n', j D = 
3d, X = 45°. 

Am'(x > 0) = fign'4- (Q — ~?'^ V exp(-x/l'), (24) 
j^, i^ + cos^x 

where jo = enl/r, I = \J Dt is spin diffusion length, v = {j' D / ]d) (N/N')\ 
the current flows along positive direction of x axis. With relationship a± = 
e 2 D±g± |14| taking into account, the parameter v can be represented as 

z = i^cFY (25) 

Quantity Z has dimensionality of contact resistance (Ohm x cm 2 ), so that the 
parameter v that determines spin current matching between two ferromagnets 
may be treated as a ratio of "spin resistances" . 

At v 3> 1, the cathode layer works as an ideal injector with equilibrium 
spin polarization (Am = 0), while equilibrium breaks in the anode layer. In 
the opposite case, v <C 1, ideal collector regime takes place, in which spin 
equilibrium breaks in the cathode layer. The nonequilibrium magnetization 
distribution at difference values of v parameter is shown in Fig. ^ 
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5 Spin accumulation resistance 



The spin equilibrium disturbance contributes to the resistance of the system in 
study. We have 

dx e dx 



J± = -<r± 

With Eqs. JSJ and (|16() taking into account, 



(26) 



dtp 
dx 



3 + 



1 



d(+ 

dx 



d(- 

dx 



e , „ „ . dm 
j + — (D + - D_) — 
2/j.b dx 



(27) 

By integrating Eq. (|27[1 over x with potential drop at x = taking into 
account, we obtain 



^(-0)-^(-L) + ^(L')-^(+0) 
3 T ' 6 -(D' + -D'_) Am' (+0), 



7 e 
-i- ~ 

(7 



cr 2/isf' 



(D+ - D_) Am(-O) 
(28) 



where L, L' are thicknesses of the contacting layers (L 3> i, £' ^ £')• 
The total potential drop over the whole system is 

U = tp(—L) - <p(L') = (^ + ^jj+ [^(-o) - vj(+0)] 

- (£>+ - D_) Am(-O) - — (£>'- Am'f+0) 



+ 



(29) 



By equating the content of the first curly brackets in Eq. I|14|l to zero, we 
find the part of the potential drop at the interface: 

p(-q) - p(+o) = (c; + c - c+ - c-) 

= ^(c' -0+4 (ac; + Ad - ac+ - ac-) = he - o 



Vm(-O) 



By substituting Eq. JSDJl into we get after some manipulations 



U = U + 



1-lB 



Am(-O) - -2— — Am'(+0) 



Z 



where 



is the voltage in absence of the spin equilibrium breaking. 



(30) 

(31) 
(32) 
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Substitution of Eqs. JSJ and into Eq. £HJ gi ves 



U = U +jZ ; o ■ ( 33 ) 

v + cos z x 



The contribution of spin accumulation to the resistance as a function of x angle 



A/?(X)/AR(JI) 




30 60 90 120 150 180 



Figure 2: Spin accumulation resistance as a function of the angle \ between 
the layer magnetization vectors at different values of the spin resistance ratio v: 
1 — v = 10 (red), 2 — v=l (green), 3 — v = 0.1 (blue). Q = 0.35, Q' = 0.15. 



ARU) = ^ = ZZ' Z^*!. . 04) 



can be found: 

u -Uo = zz , {Q^x~Q') 2 

j Z + Z' cos 2 x 

The angular dependence of AR is shown in Fig. [21 Note that AR vanishes 
at x = arccos(Q'/Q). 

The results obtained may be considered as a generalization of those in 
Refs. El IS| to the case of nonidentical noncollinear fcrromagncts. Really, 
in Ref. jH] the contact potential drop was calculated, that is, in essence, the first 
square bracket in Eq. On the other hand, in Ref. JS] the volume potential 

difference, that is the last square bracket in Eq. was calculated. As it is 
seen from Eq. I|29|l . we should sum the brackets to obtain the final result. In 
addition, we take an arbitrary angle x instead of collinear orientation taken in 
Refs. and Q3J. 
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Giant magnetoresistance (GMR) can be found from Eq. 

_ AR(n) - AR(0) 4QQ' 

= ARM " (QTW • (35) 

At Q = Q' we have Ai?(0) = 0, so that GMR takes its maximum value 
GMR = 1 . 



6 Conclusion 

We show the longitudinal spin flux continuity at the junction interface follows 
directly from the spin transformation properties under the rotation of a quan- 
tization axis. 

We derive for the first time the continuity conditions of mobile electron chem- 
ical potentials at the interface of two ferromagnetic junction layers having an 
arbitrary angle between their magnetization vectors. When the conditions were 
derived, the only statement was employed significantly, namely, the conservation 
low of the mobile electron energy flux density at the interface. 

Electron magnetization distribution in the junction is calculated based on the 
boundary conditions derived. Matching parameter is discussed, which determine 
the spin flux penetration through the interface. The parameter may be treated 
as a ratio of the contacting layers spin resistances. 

We show the disturbance of spin equilibrium at the interface leads to angle 
dependent shift of a contact potential drop and to spin accumulation magne- 
toresistance. The last effect was numerically estimated and the angles are found 
of minimal and maximal magnetoresistance. 
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